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ABSTRACT 

We present mid infrared (Mid-IR) spectra of the Compton-thick Seyfert 2 galaxy NGC3281, 
obtained with the Thermal-Region Camera Spectrograph (T-RcCS) at the Gemini South tele- 
scope. The spectra present a very deep silicate absorption at 9.7 /mi, and [Siv] 10.5 /mi and 
[Nell] 12.7 /tm ionic lines, but no evidence of PAH emission. We find that the nuclear optical 
extinction is in the range 24 < Ay < 83 mag. A temperature T = 300 K was found for the 
black-body dust continuum component of the unresolved 65 pc nucleus and at 130 pc SE, while 
the region at 130 pc reveals a colder temperature (200 K). We describe the nuclear spectrum of 
NGC 3281 using a clumpy torus model that suggests that the nucleus of this galaxy hosts a dusty 
toroidal structure. According to this model, the ratio between the inner and outer radius of the 
torus in NGC 3281 is Ro/Rd = 20, with 14 clouds in the equatorial radius with optical depth 
of Ty = 40 mag. We would be looking in the direction of the torus equatorial radius (i = 60°), 
which has outer radius of Ro ~ 11 pc. The column density is N# « 1.2 x 10 24 cmT 2 and iron Ka 
equivalent width (w 0.5 - 1.2 keV) are used to check the torus geometry. Our findings indicate 
that the X-ray absorbing column density, which classifies NGC 3281 as a Compton-thick source, 
may also be responsible for the absorption at 9.7 /mi providing strong evidence that the silicate 
dust responsible for this absorption can be located in the AGN torus. 

Subject headings: galaxies: Seyfert - infrared: ISM - ISM: molecules - techniques: spectroscopic 



1. Introduction 

Mid-infrared (Mid-IR) spectra of active galac- 
tic nuclei (AGNs) are rich in emission fea- 
tures, such as polycyclic aromatic hydrocar- 



bons (PAHs), molecular hy drogen and promi- 
ncnt forbidden emission lines (ISturm et al. Il2000t 



Weedman et al. ll2005HWu et al. ll2009HGallimore et al 
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2010l : ISales et al. 1 12010[ ) . Another spectral fea- 
ture commonly observed in AGNs is silicate at ~ 
9.7 /mi and 18 /xm, both in emission and absorp- 
tion. 

The AGN unified model proposes the existence 
of a dense concentration of absorbing material sur- 
rounding the central engine in a toroidal distribu- 
tion, which blocks the broad line region (BLR) 
emission from the line of sight in Seyfert 2 (Sy 2) 
galaxies (jAntonucci 1119931 ) . However, it is not yet 
clear if the molecular gas and dust detected in 
Mid-IR spectra is actually associated with the ab- 
sorbing material required by the so-called torus in 
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the unified model. 



HaoetaL 



(|2007l ) studied the distribution 
of the silicate feature strengths in a sample of 
196 AGNs and ultraluminous infrared galaxies 
(ULIRGs). They found that the spectra of quasars 
and Seyfert 1 (Syl) galaxies are characterized by 
silicate features in emission, with few Sy Is pre- 
senting weak absorptions. In contrast, Sy2 spec- 
tra are dominated by weak silicate absorption. 
These results suggest that the silicate emission 
(or absorption) is related to the line of sight, in 
the framework of the AGN unified model. In addi- 
tion, most of the ULIRGs sh ow very deep silicate 



absorption (|Hao et al. 1120071 ). 



The silicate behavior in AGNs can be explained 
using Nenko va's torus models (INenkova et al. 



2002 



which consider an ensemble of 
individual clouds instead of the standard "solid 
doughnut". Mason et al. ( 20061 ) detected the sili- 
cate emission feature in th e Sy 2 galaxy NGC 2110 , 



(20081)) 



and compared it with the INenkova et al. 
clumpy torus models, demonstrating that the 
presence of silicate emission can also be explained 
by an edge-on clumpy distribution. From their 
best-fit model, they conclude that there would be 
a small number of clouds along the line of sight, 
with a distribution that does not extend far above 
the torus equatorial plane. 
According 



to 



Hao et al. I (|2007l) the pres 



ence of strong silicate absorption and the ab- 
sence of PAH features in some active galaxies, 
can be interpreted as extremely heavily obscu- 
ration in these sources. An example of such 
an object is NGC 3281. Its spectrum shows a 
very deep silicate absorption at 9.7 /mi, as well 
as the forbidden emission lines [S in] 18.7 /im, 
[S iv] 10.5 /mi, [Nc n] 12.8 /im, [Nc in] 15.5 /im, 
[Ncv] 14.3/mi and 24.3 /mi, [Ncvi]7.6/mi and 
[O iv] 25.8 /im. No PAH emission was detected 



in the Spitzer spectrum ( Wu et al~l 20091 ) . In ad- 



dition, NGC 3281 is a luminous infrared galaxy 
(logL/fl ~ 10.73 L ? . ISanders et al.l 120031 ). classi- 



fied as a Sy2. Its optical spectr um was studied by 
IStorchi-Bergmann et al. ( 1992 ). finding that the 
kinematics of the kpc extended ionized gas can 
be described by rotation in a plane, with mate- 
rial outflowing (v ~ 150 km s -1 ) from the nucleus 
within a conical structure. They also suggest that 
NGC 3281 has dust asymmetrically distributed 
with respect to the ionized gas. 



In addition, using Advanced Satellit e for Cos- 
molog y and Astrophysics (ASCA) data. Isimpson I 
(|l998l ) found that NGC 3281 shows an extremely 
complex X-ray spectrum, with a prominent iron 
emission line and a large absorbing column den- 
sity (Nh ~ 7.1 ± 1.2 x 10 23 cm" 2 ) that when com- 
bined with the nuclear extinction inferred from the 
near-infrared color analysis (Ay = 22 ±11 mag), 
results in a N#/Ay ratio one order of magnitude 
highe r than that of the Galaxy ( Bohlin et al. I 
Il978l ). This effect was attributed to an opti- 
cally thick material along the line of sight of 
both the X-ray and the infrared emitting re- 
gions, which obscures the entire X-ray source but 
only a fraction of the much larger IR emi tting 
region. Later, IVignali fe Comastri ] (120021) us- 
ing the Italian-Dutch satellite BeppoSAX, found 
a Nff /Ay ratio about twice that derived by 



Simpson I (|19981 ). adding NGC 3281 to the list 



of AGNs presenting dust properties different 
from the Galaxy. The difference in the N#/Ay 
values derived by both studies, according to 



Vignali fc Comastri I (120021) . is due to the l imited 
ASCA band-pass. IVignali fc Comastri I (|2002f ) 



also show that its X-ray spectrum (E<10keV) 
is reflection-dominated, with a relatively strong 
Kq iron emission line and a heavily absorbed nu- 
clear continuum (N# ~ 1.5 — 2 x 10 24 cm~ 2 ). 
Hence, they classify NGC 3281 as a Compton- 
thick galaxy. The high values found for the ab- 
sorption column density in NGC 3281 also ex- 
plain the extremely high ratio ^[oiii] /^2-wkeV 
0.3, when compared to the average 0.02) 
value obtained for a sample of Sy2 galaxies 
(iMulchaev et al. Ill994l ). 



It follows from the above discussion that this 
galaxy has a very heavily obscured nucleus, which 
makes it a key object to study whether or not 
the observed dust absortion is associated with 
the dusty torus of the unified model. Thus, it 
may provide unambiguous observational evidence 
of this structure. We present in this paper ground 
bases Mid-IR high-angular resolution spectra of 
the galaxy NGC 3281. The Thermal-Region Cam - 



era Spectrograph (T-ReCS: iTelesco et al. Ill998[ ) 



attached to the Gemini South telescope give a spa- 
tial resolution of 16 pc/pix, which is very adequate 
to study the dust distribution in the ~ 200 pc cen- 
tral radius of this galaxy. This paper is structured 
as follows: in Section [2] we describe the observa- 
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tion and data reduction processes. Results are dis- 
cussed in Section [3] Summary and conclusion are 
presented in Section |U 

2. Observations and Data Reduction 

Ground-based Mid-IR observations of NGC 3281 
were obtained with the T-ReCS in queue mode at 
Gemini South, in 2009 April 06 UT, as part of 
program GS-2009A-Q-34. Conditions were clear, 
with water vapor column in the range 5-8mm, 
and image quality of 0'.'33 in the N band, as mea- 
sured from the acquisition images of the telluric 
standards observed right before and/or after the 
science frames. All observations used a standard 
chop-nod technique to remove time-variable sky 
background, telescope thermal emission, and the 
effect of 1/f noise from the array/electronics. The 
slit was oriented along PA = 315°, with a chop 
throw of 15", oriented along PA = 225°, perpen- 
dicular to the slit direction, in order to include 
only the signal of the guided beam position in 
the frame and avoid possible nod effects in the 
spatial direction. The same slit position/nod ori- 
entation were used for the telluric standards. The 
instrument configuration used the low-resolution 
(R ~ 100) grating and the 0'.'36 slit, for a dis- 
persion of 0.0223 /mi/pixel and a spectral reso- 
lution of 0.08 fim. The pixel scale is 0.089 arc- 
sec/pixel in the spatial direction, and the slit is 
21'.'6 long. Spectral coverage was of 7.1/im cen- 
tered at 10.5/im. The total on-source integration 
time was 980 sec. 

The data reduction was performed using the 
MIDIR and GNIRS sub-packages of Gemini irafQ 
package. To extract the final spectrum, we com- 
bined the chop- and nod-subtracted frames using 
the tasks tprepare and mistack in the midir 
package. Wavelength calibration was obtained 
from the sky lines in the raw frames. To remove 
the telluric absorption lines from the galaxy spec- 
trum, we divided the extracted spectrum for each 
observing night by that of the telluric standar d 
stars HD3438 or HD4786 (|Cohen et al. Ill999l) . 



observed before and/or after the science target. 
Finally, the spectrum was flux calibrated by inter- 



X IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Uni- 
versities for Research in Astronomy (AURA) under coop- 
erative agreement with the National Science Foundation. 



polating a black-body function to the spectrum of 
the telluric standard. For this step we used the 
task MSTELLURIC in the MIDIR package. 

Fig. [T]a, shows the slit superposed on the 
NGC 3281 N-band acquisition image. Panel b of 
this figure illustr ates the slit on the nar row [O ill] 
maps taken from lSchmitt et al.l (|2003bl ). Panel c 
shows that the N-band galaxy luminosity profile 
is extended with respect to that of the N-band 
stellar profile of HR3438. 

We extracted seven one dimensional (ID) spec- 
tra: one centered in the unresolved nucleus (4- 
pixel = 0.36 arcsec, which corresponds to ~65 
pc for a distance of 43 Mpc, using H = 74 km 
s _1 Mpc -1 ), plus extractions centered at 130 pc, 
195 pc and 260 pc to the Northwest (NW) and 
Southeast (SE) directions. Deep silicate absorp- 
tion is observed at the nucleus and at 130 pc in 
both NW and SE directions of the nucleus (Fig. 
[2]). The extraction at 195 pc NW shows a weak 
absorption. The remaining extraction does not 
present silicate absorption. Therefore, we can con- 
clude that the dust is concentrated inside a radius 
of 200 pc. 

3. Results and Discussion 
3.1. The Mid-IR Spectra 

The NGC 3281 T-ReCS spectra (Fig. clearly 
show that a deep silicate absorption, as well as the 
[Siv]10.5/xm and [Ne II] 12.7 /xm emission lines, 
are observed up to 130 pc from the nucleus. In ad- 
dition, PAH emission bands are not detected in the 
spectra, in agree ment with the Spitzer spectrum 



(|Wu et al. I [20091 ). As discussed in the introduc- 



tion, NGC 3281 is one exception among Sy 2 galax- 
ies, as most (~80%) of the objects in this class 
present PAH emission bands and no deep silicate 



absor ption in their Mid-IR spectra (pales et al. 
20101) . 



In order to estimate the silicate absorption ap- 
parent strength (S s u), we use the definition of 
Spoon et al.l(|2007l) : 



f obs (9.7 fim) 
fcont (9.7 fim) ' 



(1) 



where fobs is the observed flux density and fcont 
is the observed mid-line pseudo-continuum flux 
density. To determine the S s u, we avoid telluric 
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bands and emis sion line regions, us ing the same 



methodology as iMason et al. I (|2006l ). The values 



derived from the four spatial positions where the 
silicate absorption is observed are shown in col- 
umn 7 of Tab. [Q 

For the unresolved nucleus, we obtain S s u = 
— 1.3 ± 0.1. Similar values of S s j; are observed in 
ULIRGs, while much weaker (average i s -0.61) ab 



sorpt ions are seen in Sy 2s (see Fig. 2 of lHao et al. 



20071 ). The remaining three extraction have simi- 



lar values to that of the nucleus. 

When dealing with dusty sources, a more phys- 
ical parameter is the line of sight depth (r), 
which can be estimated from S s u at 9.7 ^m us- 
i ng the equation ln(e~ T 9 7 ) = -r 9 . 7 = 
( Nenkova et al. ~l l2008bl) . furthermore, we can de- 
rive the apparent o ptical extinctio n using A^P = 
T9.7 x 18.5 ±2 mag (jPraine Il200lh . The resulting 
values are shown in Tab. Q] The dust extinc- 
tion for the unresolved nucleus of NGC3281 is 
papp = 24 ± 5 mag, in agreement with the value 
(Ay = 22 ± 11 m ag) derived by combined IR data 
( Simpson 19981) . Peruse that our uncertainties is 
much lower than these authors. We note from the 
values listed in Tab. [T] that the position labeled 
as 1 (130 pc from the nucleus at NW direction) 
presents a slightly higher extinction than the cen- 
tral one (Ay = 28 ±8 mag). 

To estimate in a more robust way the line 
of sight depth at 9.7it m (79.7) in NGC 3 281 we 
used the PAHFItH code ( Smith et al. 112007 ). which 
assumes that the source spectrum is composed 
by continuum emission from dust and starlight, 
prominent emission lines, individual and blended 
PAH emission bands, and that the light is atten- 
uated by extinction due to silicate grains, pah- 
fit m odels the extinc t ion us ing the dust opacity 
law of iKemper et al. I ([20041 ) . where the infrared 
extinction is considered as a power law plus sili- 
cate features peaki ng at 9.7 /zm. For details see 
Smith et al. I (j2007l ). 



We used this approach to determine the con- 
tribution of the distinct components to the spec- 
tral energy distribution (SED) of NGC 3281 as a 
function of the distance to the galaxy centre. As 
the spectra do not show PAH emission, we have 
not included this component in the fitting. The 



ot her input par a meter s are the same as those used 
bv lSales et al~l (|2010l ). which arc appropriate for 
AGNs. The results of the decomposition of the 
NGC 3281 spectra are shown in Fig. and the 
derived silicate optical depths are given in Tab. [1] 

Clearly, T9.7 varies along the slit, being T9.7 
= 4.5±0.7 at the unresolved nucleus and t 9 7 = 
4.6±0.7 at 130 pc NW extraction. In constrast, it 
reaches T9.7 = 5.5±0.8 for the 130 pc SE extrac- 
tion. The dust extinction (Ay) values estimated 
using the PAHFIT derived line of sight depth, col- 
umn 10 in Tab. [TJ are much larger than those ob- 
tained from the S s u . The latter are obtained from 
the depth at 9.7/im while the pahfit code takes 
the whole profile of the silicate absorption fea- 
ture into account, thus the S s j; indicator should 
be taken as a lower limit for the Ay. The ob- 
tained opacity values we re compared to those o f 
the Seyfert sample from iGallimore et al. ( 2010h . 
The histogram of Tg.7 (Fig. [5]) shows that the dust 
opacity in NGC 3281 is one of the highest observed 
in Seyfert galaxies. 

The thermal dust continuum components re- 
quired to fit the observed spectra of the unresolved 
nucleus and 130 pc SE is a black-body with T = 
300 K, while the 130 pc NW spectrum requires it 
a black-body of T = 200 K. This colder region 
has the deepest silicate absorption, consequently is 
more heavily obscured (Ay = 102 ± 26 mag) than 
the nucleus and the SE region. The presence of 
this obscuring material is consistent with the spa- 
tial distribution of the optical reddening found by 



Storchi-Bergmann et al.l (|1992l ) . Note that our slit 



position is perpendicular to the cone observed in 
[Om]. See their Fig. 9 and 2. 

In order to get additional support of that highly 
reddened region, we overplot the T-RcCS slit 
position on t h e [Oin lA5007A maps taken from 
Schmitt et all (|2003bh . We found that the 130 pc 



NW region coincide with a region where [Om] 
emission is not observed in agreement with those 
previously found from X-ray data. That is, the 
NGC 3281 is a heavily obscured source. 



2 Source code and documentation for PAHFIT are available at 
http: / /tir. astro. utoledo.edu/jdsmitli/rcscarcli/palifit.plip 
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3.2. Is the observed silicate absorption de- 
tected at the unresolved nucleus asso- 
ciated with the dusty torus required 
by the unified AGN model? 

There are several torus models, some authors 
take it as h aving an uniform dust density distri 
bution (e.g. Pier fc Krolik Il992t Granato et al. 



1997[ ISiebenmorgen et al. I 12004 iFritz et al. 
2006). However, for dust grains to survive in 



the torus environmen t, they should be forme d 
by clumpy structures ( Krolik fc Begelmanl 1988 ). 
and new models have been developed where the 
dust grains are distributed in a c lu mpy medium 
(e.g jNenkova et al. I [20021 l2008allbl ; iHonig et al." 
20061 ). Therefore, we compare the uncontami 



nated nuclear spectrum (see below) with theoret- 
ical SEP obtained f r om Ne nkova's models (e.g. 



Nenkova et al. 1120021 12008c 

In their modeling, Nenkova et al. assume that 
the dusty material surrounding the AGN is dis- 
tributed in clumps forming a toroidal geometry, 
and determine the following parameters: i) the 
number of clouds Nq in the torus equatorial ra- 
dius; ii) the optical depth for an individual cloud 
in the V band (jv)] Hi) the radial extension of the 
clumpy distribution, Y = Ro/Rd, where Rq and 
Rd are the outer and inner radius of the torus, 
respectively; iv) the power law index q that de- 



scribes the radial density profile (oc 



v) the 



torus angular width, constrained by a Gaussian 
angular distribution of cr; vi) the observed view- 
ing angle i. 

In order to properly analyze the spectrum from 
the nuclear dusty structure, it is necessary to first 
remove the host galaxy contribution. For this, 
we have subtracted from the central aperture the 
spectrum obtained from averaging the two adja- 
cent extractions (apertures -1 and 1, Fig.[T]). Noise 
effects arc reduced by smoothing (with a rectan- 
gular box of size 5 pixels) the final spectrum. The 
emission lines were removed by simple interpola- 
tion and the region where the telluric band is lo- 
cated (Figj2) were not used in the fit. Then, the 
smoothed spectrum was compared with the ~ 10 6 
clumpy theoretical SEDs, and the best fit is ob- 
tained by searching for the minimum in the equa- 
tion: 



3 The models are available 

http: / /www.pa.uky.cdu/clumpy/ 



from 



1 N 

-Y 

i=l 



F 



obs. Xi 



F 



mod, Xi 



'A, 



(2) 



where N is the number of data points in the spec- 
trum, F f, Sj Xi ' Fmod. Xi and 5\ i are the observed, 
theoretical fluxes at each wavelength and their re- 
spective uncertainties. The later were taken as 
being 10% of Fa, which are the expected ones 



for T-ReCS observ ations (jRadomski et al. 112002 : 



Mason et al. 1 120061 ) . Note that both F bs,\i and 
Fmod, \i were normalized to unit at 9.0/zm, with 
the uncertainties being properly propagated. The 
parameter set that provides the best fit is shown 
in Table [2] and the corresponding theoretical SED 
is over-plotted on the NGC3281 decontaminated 
nuclear spectrum (Fig. 2]). Due to the very large 
number of models required to cover the parameter 
space, we estimate the uncertaintie s for the bes t 
model followi ng the approach of iMason et al. 



(|2009i see also lNikutta et al. I (|2009h ). but instead 
of using only the three best models, we calculated 
the mean and standard deviation for the parame- 
ters considering all models with a \ 2 value within 
10% of the best-fitting result (76 in total). The lo- 
cus of the corresponding theoretical SEDs is plot- 
ted as a grey region in Fi gure IH Our approac h 



is similar to that used by iNikutta et al. 



thus, we have also tested other y 2 deviation frac- 
tions (5%, 10%, 15% and 20%, corresponding to 
17, 76, 151 and 291 acceptable solutions, respec- 
tively). In general our technique produce simi 



l ar re sults than that obtained by INikutta et al. 
(|2009l ). rowever, our methodology tends to pro- 



duce a locus with less acceptable solutions than 
Nikutta's one. The mean parameters derived with 
the different \ 2 deviations are shown in Table [2] 

The physical parameters of the best-fit model 
to the nuclear uncontaminatcd SEDs suggest that 
NGC3281 hosts a dusty toroidal structure. The 
dusty clouds, each with an optical depth ry = 
40 mag, occupy a toroidal volume with Rq /Rd = 
20. The distribution of the clouds follows a power- 
law radial behavior (r -15 ), and the number of 
clouds along the equatorial radius is 10. The an- 
gular distribution of the clouds is characterized by 
a width a = 70°. These results are in agreement 
with the unified AGN model, which requires the 
presence of a dusty structure that obscures the 
broad line region in Sy 2 galaxies. From the fitted 
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models we would be looking in a direction nearby 
of the torus equator (i = 60°) 

In general, the torus para meters derived here 
are sim ilar to those obtained bv lRamos Almeida et al. 
( 20091 ), except for the number of clouds and op- 
tical depth per single cloud, where they find only 
5 clouds with ry = 10 in the equatorial radius. 
This discrepancy can be due to the different 
constraints available, 



as 



Ramos Almeida et al. 



(2009 ) used Mid-IR photometric data (see also 



Ramos Almeida et al. 1201 it lAlonso-Herrero et al. 



20111 ). while the data presented here would allow 
for a more detailed description of the spectral be- 
havior of the silicate dust, thus providing a better 
input to the modeling procedure. 



The Nenkova et al. 



(|2008bl ) models adopt a 
standard dust-to-gas ratio, and assume the col- 
umn density of a single cloud to be ~ 10 22 — 
10 23 cm~ 2 . In addition, the total number of clouds 
at a viewing angle i is given by a Gaussian distri- 



bution: N b s (i) — N exp 



thus result- 



ing in a column density of N^ bs ' > = N b s (i)NH 
at the observer view angle. Using the above, 
we can then derive a column density of N# w 
1.2 x 10 24 cm~ 2 for the uncontaminated nuclear 
spectrum of NGC 32 81, which is consistent with 
the values derived bv IVignali fc Comastri I ( 2002 ) 
from X-ray observations (~ 2 x 10 24 cm -2 ). 



Moreover jLevenson et al.~ ( 2002 ) demonstrated 
that for heavily obscured AGN the iron Kq 
EW values depend on the torus geometry. We 
compare the torus geometrical parameters de- 



rived in our paper with those of iLevenson et al 
(|2002h . and found that the present torus aper- 
ture (a = 70°) corres pond to = 20° in the 



Levenson et al. I (|2002l )'s models. The observer 
view angle has the same definiti on for both mod 
els_(i = 60°). From Figure 2 of ILevenson et al. 



((2002J) and using the above parameters we esti- 
mate that the iron Ka EW for NGC 3281 is w 
2 - 3keV, which is in ag r eeme nt with that found 



by IVignali fc Comastri ( 2002 ). confirming their 



results, that this galaxy would shows a large EW 
for the iron Ka line. We point out that this torus 
geometry was inferred from the silicate absorp- 
tion line, this indicates that the X-ray absorb- 
ing column density, which makes NGC 3281 a 
Compton-thick source, may also be responsible 
for the absorption at 9.7 /iin. Therefore, our re- 



sults provide strong evidence that the silicate dust 
responsible for such abs orption is located in th e 



AGN torus. Following iNenkova et al. I (|2008bl ). 



and adopting a dust sublimation temperature of 
1500K, we estimate a torus outer radius of Ro ~ 
11 pc. This value is consistent with the results 



of .Taffe et al. (2004) 



Ramos Almeida et al. 



Tristram et al. 
(|2009l) . 



(120071 ) and 



We also derive the intrinsic bolometric lu- 
minosity for the AGN from the best-fit model, 
finding L bo ; = 1.9xl0 45 erg s -1 . This value 

ag rees with the intrin sic luminosity obtained 

bv IVignali fc Comastri I (jiooi) from X-ray data, 
when applyi n g the conversion factor suggested by 



lElvis et al. 



( 20041 ) . resulting L b °'_ ray « 3.2 x 
10 44 erg s" 1 . We note that NGC 3281 has a high 
X-ray luminosity when compa r ed to objects stud- 



ied by Ramos Almeida et al. ( 20091 ). 



4. SUMMARY AND CONCLUSIONS 

In this work we present high-spatial resolution 
Mid-IR spectra of the Compton-thick Sy 2 galaxy 
NGC 3281. The ground-based Mid-IR data were 
taken with the T-ReCS spectrograph, and our 
main results can be summarized as follows. 

• NGC3281's spectra show a very deep sil- 
icate absorption, as well as [Siv]10.5^m 
and [Neil] 12.7/iim ionic lines. However, the 
spectra do not present PAH emission bands, 
which makes t his galaxy uncomm on among 



the Sy 2s class (|Sales et al. Il2010l) . Further 



more, the dust is concentrated inside a ra- 
dius of 200 pc. 

• The optical extinction, Ay, of the nucleus 
and its neighborhood was estimated using 
the S s u indicator and pahfit code. From 
the Ssu indicator we infer the nuclear dust 
extinction Ay PP = 24 ± 2 mag, while the 
pahfit code provides Ay = 83 ±22 mag. 
We point out that the first value can only 
be used as a lower limit. This large extinc- 
tion confirms that th i s gala xy, if compared 
with lGallimore et al.l ( 2010l )' Seyfert sample 
(83 sources), is a heavily obscured source 
compatible with a Compton-thick galaxy 
scenario. 

• The black-body temperature due to dust 
continuum components for the unresolved 
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nucleus and the region at 130 pc SE is T = 
300 K. A colder black-body (T = 200 K) was 
found for the region at 130 pc NW, which 
also has the highest derived visual extinction 
(Ay = 102 ±0.2 mag). 

The nuclear uncontaminatcd spectrum fol- 



lows a clumpy torus model (INcnkova et al 



2008bh . which suggests that NGC 3281 has a 
dusty toroidal structure. The dusty clouds 
occupy a toroidal volume with Rq/R^ = 20, 
with each cloud having an optical depth ry 
= 40 mag. The distribution of the clouds fol- 



lows a power-law radial behavior 



-1.5 



and 



the number of clouds in the equatorial radius 
is 14. The angular distribution of the clouds 
has width of a = 70°. These results agree 
with the unified AGN model that to block 
the broad line region emission of Sy 2 galax- 
ies, requires a dusty structure. According to 
the best-fit model, we would be looking in 
the direction nearby of the torus equatorial 
radius (i = 60°), having an outer radius of 
Ro ~ 11 pc. 

We derive a column density of Njj = 
1.2 x 10 24 cm" 2 for NGC 3281 indicates that 
the X-ray absorbing column density, which 
classifies NGC 3281 as a Compton-thick 
source, may also be responsible for the ab- 
sorption at 9.7 fim. In addition, our results 
provide strong evidence that the silicate dust 
responsible for such absorption is located in 
the AGN torus. From Nenkova's models 
we find that the bolomctric luminosity of 
this AGN is L boi = 1.9 xl0 45 crg s _1 , which 
shows that NGC 3281 is a very luminous 
source. 
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Fig. 1. — (a) The T-ReCS long-slit position overplotted on the contoured NGC3281 acquisition image. The 
N-band emission peak is centered in the image, which was supposed to coincide with the active nucleus. 
The c ontours are linear an d stepped of 10% of the peak, (b) the slit superimposed on a [Oiii]A5007A 
maps (jSchmitt et al.ll2003bh and spectral extraction positions were labeled, (c) the T-ReCS N-band spatial 
emission profile, solid and dotted lines represent the NGC3281 galaxy and the HR3438 standard star, 
respectively. The fluxes are normalized to the peak value. 
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Fig. 2. — Spectra of NGC 3281 extracted in 65 pc steps along P.A. = 315°. Extraction positions are labeled, 
and dashed lines indicate the positions of the [Siv] 10.5 fim and [Nell] 12.7/zm ionic lines. Telluric 3 band 
is represented by the hatched area. ^ 




Fig. 3. — (a) The PAHFIT decomposition spectra of the 130 pc SE (-1) extraction, (b) and (c) are the same 
for the unresolved nucleus (0) and 130 pc NW (1) extractions, respectively. Measurements are represented 
by squares, uncertainties are plotted as vertical error-bars, which are smaller than the symbol size. Dotted 
black line indicates mixed extinction components. Solid grey and black line represent total and individual 
thermal dust continuum components, respectively. Green lines are emission lines. Magenta line represents 
the best fit model. 
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8.5 9.0 9.5 10.0 10.5 11.0 11.5 12.0 



Fig. 4. — Red and grey lines are the best and 10% of best-fitting SEDs to the NGC3281 uncontaminated 
spectrum, respectively. All spectra were normalize in flux at 9.0/im. 



Table 1 
Line Measurements 



Position 


Label 


[S IV] 10.5 M m 


EW[Siv] 


[Nell] 12.7 fim 


EW[Neii] 


S s u 


AS?" (mag) 


TQ.7 


Av(mag) 


130 pc SE 


-1 


4.3±0.7 


0.7 






-1.4±0.1 


26 ±5 


5.5±0.8 


102 ± 26 


Center 





6.2 ± 1.0 


0.1 


1.3±0.3 


0.02 


-1.3±0.1 


24 ±5 


4.5±0.7 


83 ±22 


130 pc NW 


1 


2.3±0.4 


0.2 


2.2±0.4 


0.06 


-1.5±0.3 


28 ±8 


4.6±0.7 


85 ±22 


195 pc NW 


2 










-1.1 ±0.1 


21 ±4 







NOTE. — Columns 3 and 5 are values of the emission line integrated fluxes given in units of 10 — 16 W m~ 2 . Columns 4 and 6 are values 
of the equivalent widths in units of fim. Column 7 is the silicate strengths computed from cq. Column 8 shows the apparent optical 
extinction values derived from S s ii values. Column 9 is the values of optical depth inferred by PAHFIT and column 10 their optical 
extinction. In order to determine the optical extinction values, we assume A^ pp / r s n — 18.5 ± 2 (Drainc 2003). 
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g 35 

cr 

<U 30 - 



n 



3-4 4-5 6-7 

Dust Opacities 



Fig. 5. — Frequency histogram of the dust opacity of the 83 Seyfert galaxies taken from iGallimore et al 
(l2010h . Dotted line shows the value of the NGC3281. 



Table 2 

Fitted Parameters with CLUMPY Model 



Parameter 


Best Fit 


Average 


Average 


Average 


Average 






5% 


10% 


15% 


20% 


Torus angular width (cr) 


70° 


68° ±3° 


62° ±6° 


62° ±7° 


61° ±8° 


Radial extent of the torus (Y) 


20 


20 ±0 


23 ±4 


28 ±18 


38 ±27 


Number of clouds in the torus equatorial radius {No) 


14 


13 ± 1 


13 ±2 


13 ±2 


13 ±2 


Power-law index of the radial density profile (q) 


1,5 


1.5 ±0 


1.4±0.2 


1.4 ±0.3 


1.2 ±0.5 


Observers' viewing angle (i) 


60° 


69° ± 11° 


76° ± 11° 


76° ± 12° 


72° ± 18° 


Optical depth for individual cloud (tv) 


40 mag 


40 ± mag 


40 ± mag 


40 ± 3 mag 


41 ± 8 mag 


Total number of solution 




17 


76 


151 


291 



Note. — Columns 2. 
respectively. 



3, 4, 5 and 6 show the best value and the average for 5%, 10%, 15% and 20% of variation of the 



12 



Levenson, N. A., Krolik, J. H., Zycki, P. T., Heck- 
man, T. M., Weaver, K. A., Awaki, H., & 
Terashima, Y. 2002, ApJ, 573, 81L 

Levenson, N. A., Sirocky, M. M., Hao, L., Spoon, 
H. W. W., Marshall, J. A., Elitzur, M., & 
Houck, J. R. 2007, ApJ, 654, 45L 

Marshall, N., Warwick, R. S., & Pounds, K. A. 
1981, MNRAS, 194, 987 

Mason, R. E., Geballe, T. R., Packham, C, Lev- 
enson, N. A., Elitzur, M., Fisher, R. S., & Perl- 
man, E. 2006, ApJ, 640, 612 

Mason, R. E., Levenson, N. A., Shi, Y., Packham, 
C, Gorjian, V., Cleary, K., Rhee, J., & Werner, 
M. 2009, ApJ, 693, 136L 

Mulchaey, J. S., Koratkar, A., Ward, M. J., Wil- 
son, A. S., Whittle, M., Antonucci, R. R. J., 
Kinney, A. L., & Hurt, T. 1994, ApJ, 436, 586 

Nenkova, M., ivezic, Z., & Elitzur, M. 2002, ApJ, 
570, L9 

Nenkova, M., Sirocky, M. M., Ivezic, Z., & Elitzur, 
M. 2008, ApJ, 685, 147 

Nenkova, M., Sirocky, M. M., Nikutta, R., Ivezic, 
Z., & Elitzur, M. 2008, ApJ, 685, 160 

Nenkova, M., Sirocky, M. M., Nikutta, R., Ivezic, 
Z., & Elitzur, M. 2010, ApJ, 723, 1827 

Nikutta, R., Elitzur, M., & Lacy, M. 2009, ApJ, 
707, 1550 

Phillips, M. M., Charles, P. A., & Baldwin, J. A. 
1983, ApJ, 266, 485 

Pier, E. A., & Krolik, J. H. 1992, ApJ, 401, 99 

Radomski, J. T., Piha, R. K., Packham, C, Tele- 
sco, C. M., & Tadhunter, C. N. 2002, ApJ, 566, 
675 

Ramos Almeida, C, Levenson, N. A., Rodriguez 
Espinosa, J. M., ct al. 2009, ApJ, 702, 1127 

Ramos Almeida, C, Levenson, N. A., Alonso- 
Herrero, A., et al. 2011, ApJ, 731, 92 

Roche, P. F., Packham, C, Aitken, D. K., & Ma- 
son, R. E. 2007, MNRAS, 375, 99 



Sales, Dinalva A., Pastoriza, M. G., & Riffel, R. 
2010, ApJ, 725, 605 

Sandagc, A. 1978, AJ, 83, 904 

Sandage, A., & Brucato, R. 1978, AJ, 84, 472 

Sanders, D. B., Mazzarella, J. M., Kim, D. C, 
Surace, J. A., & Soifer, B. T. 2003, AJ, 126, 
1607 

Schmitt, H. R., Ulvcstad, J. S., Antonucci, R. R. 
J., & Kinney, A. L. 2003, ApJS, 132, 199 

Schmitt, H. R., Donley, J. L., Antonucci, R. R. J., 
Hutchings, J. B. & Kinney, A. L. 2003, ApJS, 
148, 327 

Siebcnmorgen, R., Kriigel, E., & Spoon, H. W. W. 
2004, A&A, 414, 123 

Simpson, C. 1998, ApJ, 509, 653 

Smith, J. D. T., Drainc B. T., Dale D. A., ct al. 
2007, ApJ, 656, 770 

Spoon, H. W. W., Marshall, J. A., Houck, J. R., 
Elitzur, M., Hao, L., Armus, L., Brandl, B. R., 
& Charmandaris, V. 2007, ApJ, 654, 49 

Storchi-Bcrgmann, T., Wilson, A. S., & Baldwin, 
J. A. 1992, ApJ, 396, 45 

Sturm, E., Lutz, D., Tran, D., Fcuchtgruber, H., 
Gcnzel, R., Kunze, D., Moorwood, A. F. M., & 
Thornlcy, M. D. 2000, A&A, 358, 481 

Tclcsco, C. M., Piha, Robert K., Hanna, Kevin 
T., Julian, Jeffrey A., Hon, David B., & Kisko, 
Thomas M. 1998, Proc. SPIE, 3354, 534 

Tristram, K. R. W., Mcisenhcimer, K., Jaffc, W., 
et al. 2007, A&A, 474, 837 

Vignali, C, & Comastri, A. 2002, A&A, 381, 834 

Young, S., Packham, C, Mason, R. E., Radomski, 
J. T. & Telesco, C. M. 2007, MNRAS, 378, 888 

Weedman, D. W., Hao, Lei, Higdon, S. J. U., et 
al. 2005, ApJ, 633, 706 

Wu, Y., Charmandaris, V., Huang, J., Spinoglio, 
L., & Tommasin, S. 2009, ApJ, 701, 658 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 



13 



